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We study a nonlinear g-voter model with stochastic noise, interpreted in the social context as 
independence, on a duplex network. To study the role of the multi-levelness we propose three 
methods of transferring the model from a mono- to a multiplex network. They take into account 
two criteria - one related to the status of independence (LOCAL vs. GLOBAL) and one related to peer 
pressure (AND vs. OR). In order to examine the influence of the presence of more than one level in the 
social network, we perform simulations on a particularly simple multiplex - a duplex clique, which 
consists of two fully overlapped complete graphs (cliques). Solving numerically the rate equation 
and simultaneously conducting Monte Carlo simulations, we provide evidence that even a simple 
rearrangement into a duplex topology may lead to significant changes in the observed behavior. 
However, qualitative changes in the phase transitions can be observed for only one of the considered 
rules - LOCAL&AND. For this rule the phase transition becomes discontinuous for q — 5, whereas for 
a monoplex such a behavior is observed for q = 6. Interestingly, only this rule admits construction 
of realistic variants of the model, in line with recent social experiments. 


I. INTRODUCTION 

Opinion dynamics is one of the most investigated sub¬ 
fields of sociophysics Subjectively, there are at 

least two important reasons why physicists study this 
topic. The first motivation comes from social sciences 
and can be described as a temptation to build a bridge 
between the micro and macro levels in describing social 
systems. Traditionally, there are two main disciplines 
that study social behavior - sociology and social psychol¬ 
ogy. Although the subject of the study is the same for 
both disciplines, the usually taken approach is very dif¬ 
ferent. Sociologists study social systems from the level 
of the social group whereas social psychologists concen¬ 
trate on the level of the individual [5|- From the physi¬ 
cist’s point of view this is similar to the relationship 
between thermodynamics and statistical physics. This 
analogy raises the challenge to describe and understand 
the collective behavior of social systems (sociology) from 
the level of interpersonal interactions (social psychology). 
The second motivation to deal with opinion dynamics is 
related to the development of non-equilibrium statisti¬ 
cal physics. Models of opinion dynamics are often very 
interesting from the theoretical point of view [6]. One 
of the best examples is the famous voter model or, the 
recently introduced, nonlinear q- voter model [?J. Both 
models are based on dichotomous opinions and belong to 
a wide class of binary-state dynamics mm- It should 
be stressed here, that binary opinions are natural from 
the social point of view, since dichotomous response for¬ 
mat with 1 (yes, true, agree) and 0 (no, false, disagree) 
as response options is one of the most common in social 
experiments [IjJ [l2j . 

Among many other binary opinion models M , the q- 
voter model is not only interesting from theoretical point 
of view, but also justified from the social point of view. 
In short, within the q -voter model each individual in¬ 
teracts with a set of q neighbors (a g-lobby) and if all 
q neighbors share the same state (i.e. the g-lobby is 


unanimous), the individual conforms to this state. As 
originally proposed, in the other case (disagreement) the 
individual changes its state with probability e [7]. How¬ 
ever, in some later publications the model with e = 0 was 
studied, as a natural generalization of the Sznajd model 
pMl4l5j. The unanimity rule can be justified based on 
social experiments. It has been observed in number of 
experiments that a small unanimous group may be more 
efficient than a much larger group with a non-unanimous 
majority [E[. In a classical series of experiments on con¬ 
formity, Solomon Asch has found that the presence of a 
social supporter reduced conformity dramatically - par¬ 
ticipants of the experiment were far more independent 
when they were opposed by a seven person majority and 
had a partner sharing the same opinion than when they 
were opposed by a three-person majority and did not 
have a partner [16|. Influence of a consistent minority 
on the responses of a majority has been reported by 
Moscovici et al. 0 - Also recent neurological experi¬ 
ments suggest that unanimous opinions may be critical 
for normative influence [l8]. 

From the physicist’s point of view the q -voter model is 
interesting because of the rich behavior related to phase 
transitions @,0,0 as well as the controversy related to 
the exit probability of the model 00- In this paper 
we will focus on phase transitions driven by stochastic 
noise which - in the social context - may be interpreted 
as independence [l3|, 20]. In social psychology, indepen¬ 
dence is recognized as one of the two types of noncon¬ 
formity and means resisting influence [2lj. It has been 
noticed that independence plays a role similar to the tem¬ 
perature and introduces order-disorder phase transitions 
[20| . Interestingly, it has been shown that in the case of 
a complete graph, the phase transition changes its type 
from continuous to discontinuous for g > 6 [13|, |20| . 

Up till now, the g-voter model has been studied on 
monoplex networks, i.e. networks that consist of only 
one level. However, as noted recently, interactions among 
individuals can be of qualitatively different nature and 
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therefore modeled by multi-level networks [22|. In the 
last two years, a lot of attention has been devoted to the 
analysis of various dynamics on multiplex networks, in¬ 
cluding diffusion process es f23l| . epidemic sprea ding [2j- 
[26] and voter dynamics |3t| . Brummitt et al. |27j have 
generalized the threshold cascade model on complex net¬ 
works [28|. In j29j they further expanded the model and 
introduced the idea of OR and AND nodes. An OR node 
is activated as soon as a sufficiently large fraction of its 
neighbors are active in at least one level. An AND node 
is activated only if in each and every layer a sufficiently 
large fraction of its neighbors are active. We will use a 
related notion of OR and AND types of social influence. 

Without doubt most real-world social networks con¬ 
sist of many levels. For instance, a student may belong 
to a network of classmates, a network of sport-club team¬ 
mates and a network of Facebook friends. The question 
is if this multi-levelness is important for the macroscopic 
(or global) properties of the social system, such as public 
opinion, or not. We try to answer this question within 
the g-voter model that takes into account two types of 
response to social influence - conformity and indepen¬ 
dence. 

The remainder of this paper is organized as follows. 
In Sec. II, we briefly recall the g-voter model on a 
single monoplex clique. We extend this framework to 
multiplex networks in Sec. Ill and propose three rules 
(GLOBAL&AND, GLOBAL&OR and LOCAL&AND). In Sec. IV 
we derive rate equations, that describe the time evolution 
of the system, for each of the three rules. Furthermore 
in Sec. V, based on these equations, we derive phase 
diagrams and compare results obtained from analytical 
equations with those obtained from Monte Carlo simula¬ 
tions. The next section is devoted to a deeper analysis 
of the phase transitions. We wrap up the results and 
conclude in Sec. VII. 


II. THE q -VOTER MODEL ON A SINGLE 
MONOPLEX CLIQUE 


The g-voter model with the stochastic noise, inter¬ 
preted in the social context as independence, has been 
already analyzed on a monoplex complete graph [l3| . 
Here, following [Tj|, we also consider a set of N individu¬ 
als, which are described by the binary variables Si = ±1 
(spins ’up’ or ’down’). At each elementary time step t we 
randomly choose an i-th node (i.e. a voter) and a g-lobby, 
which is a randomly picked group of q individuals. Only 
if the g-lobby is self-consistent it can influence the voter. 
With probability 1 — p the g-lobby (if it is homogeneous) 
acts on the state of the voter, which means that the voter 
changes state to the state of the g-lobby. With probabil¬ 
ity p voter behaves independently - with equal probabil¬ 
ities flips to the opposite direction Si(t + 1) = —Si(t) or 
keeps its original state Si(t + 1) = Si(t). Therefore only 
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where a single-line arrow represents the state of a node 
belonging to the g-lobby while the state of the voter is 
marked with a double-line arrow. It has been shown that 
the system, described by the g-voter model with indepen¬ 
dence, undergoes the phase transition at p = p c {q ). For 
p < p c the majority coexists with the minority opinion 
(ordered state) and for p > p c there is a status-quo (dis¬ 
ordered state) 0. Interestingly, it occurred that for 
q < 5 the phase transition is continuous, whereas for 
q > 5 becomes discontinuous. 

The most natural quantity that describes the macro¬ 
scopic behavior of such a system is magnetization, which 
from the social point of view represents so called public 
opinion: 


1 

= ( 2 ) 

v i =1 

Moreover, in the case of a complete graph, the magneti¬ 
zation fully describes the state of system. In this paper 
we will calculate it in two ways - by Monte Carlo sim¬ 
ulations of the microscopic system of the size N and by 
numerical solution of the equation describing the time 
evolution of the average magnetization. In the case of 
the Monte Carlo simulations we will calculate ensemble 
average of the magnetization in the stationary state: 

1 M 

<m>= ( 3 ) 

3 = 1 

where mi denotes the stationary value of the magnetiza¬ 
tion in i-th realization (sample) and i = 1,..., M. In this 
paper we average all Monte Carlo results over M = 10 3 
samples. 


III. THE q- VOTER MODEL ON A DUPLEX 
CLIQUE 

Let us start with defining a duplex clique, which is a 
particular case of a multiplex. Specific definitions of mul¬ 
tiplex networks have been introduced in j22|,[3l],HU . Such 
systems consist of distinct levels (layers) and the inter¬ 
connections between levels are only between a node and 
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FIG. 1. The topology of a duplex clique, i.e. a network that 
consists of two distinct levels (layers), each of which is repre¬ 
sented by a complete graph (i.e. a clique) of size N. Levels 
represent two different communities but are composed of ex¬ 
actly the same people - each node possesses a counterpart 
node in the second level. Interconnections between levels (de¬ 
noted by the gray dashed lines) are realized exclusively by 
connecting the node with its own counterpart (i.e. the same 
node) on the other level. Links within a single level (black 
solid lines) represent some kind of social relation (e.g. friend¬ 
ship) . 


FIG. 2. The topology of two inter-connected monoplex 
cliques. This network consists of two levels, each of which 
being represented by a complete graph of size N. Levels rep¬ 
resent two different communities and are composed by dif¬ 
ferent people. Links within a single level are equal to the 
connections between levels (black solid lines) and represent 
some kind of social relation (e.g. friendship). 
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its counterpart in the other layer (i.e. the same node). 
Here we consider a duplex clique (see Fig. [1]), i.e. a net¬ 
work that consists of two distinct levels (layers), each of 
which is represented by a complete graph (i.e. a clique) of 
size N. Levels represent two different communities (e.g. 
Facebook and school class), but are composed of exactly 
the same people - each node possesses a counterpart node 
in the second level. Such an assumption reflects the fact 
that we consider fully overlapping levels, being an idealis¬ 
tic scenario. We also assume that each node possesses the 
same state on each level, which means that the society 
consists of non-hypocritical individuals only. 

It is worth to stress the difference between a duplex 
clique and two inter-connected monoplex cliques (see Fig. 
0), where the state of a node on one level is not directly 
related to the state of a node in the second layer. In 
inter-connected monoplex cliques the inter-clique links 
fulfill the same role as intra-clique edges. Classical voter 
models were analyzed on inter-connected cliques in [3jl 
0 . 

In this paper we investigate the q -voter model on a du¬ 
plex clique. As in the case of the monoplex, we consider 
a set of N individuals described by the binary variables 
Si = ±1, which are the same on both levels. On each 
level all individuals are connected with each other and 
therefore create a duplex clique, as shown in Fig. [lj 

We consider two criteria of level dependence - one re¬ 
lated to the status of independence (GLOBAL vs. LO¬ 
CAL) and one related to peer pressure (AND vs. OR), 
(see also Tab. [J): 

1. Criteria related to the status of independence: the 
GLOBAL rule means that an agent is independent on 
both levels, but the LOCAL rule admits a situation 


TABLE I. Three versions of dynamics on the multiplex clique. 
We do not consider the LOCAL&OR rule because of the difficulty 
of such a concept related to social unreality and algorithmic 
ambiguity. 


where a person is independent in one clique but not 
in the other. 

2. Criteria related to the peer pressure: the AND dy¬ 
namics is more restrictive and a node changes its 
state only if both levels suggest changes, in the OR 
variant one level is enough to change the state of 
an individual. 

Finally we propose the following three rules: 

(i) GLOBAL&AND - global independence and the AND rule 
(see an example in Fig. [3]) 

With probability p the voter is independent and 
with 1 — p behaves like a conformist regardless of the 
level. In the case of independence the voter changes 
its state to the opposite one with probability 1/2 
(we automatically change the state of the voter on 
both levels). In the case of conformity the voter 
changes its state only when both ^-lobbies (i.e. on 
the first and on the second level) are homogeneous 
and both have the state opposite to the state of 
the voter. Therefore only the following changes are 
possible: 
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FIG. 3. The GLOBAL&AND rule: a voter is independent regard¬ 
less of the level with probability p and is subjected to the peer 
pressure with probability 1 — p only if (/-panels on both lev¬ 
els are self-consistent. In this example on level 1 the q -lobby 
(agents in circles) is homogeneous but on level 2 the lobby 
is not self-consistent. Therefore the voter will not change its 
state under the peer pressure. 



FIG. 4. A comparison between the (/-voter model on the 
monoplex clique (full symbols) and the (/-voter model with 
the GLOBAL&AND rule on the duplex clique (empty symbols) - 
the GLOBAL&AND rule leads to a trivial result, identical with the 
monoplex case for a doubled value of q. An ensemble average 
< m > of the magnetization, as a function of the stochastic 
noise p, was obtained by the Monte Carlo simulations for the 
system of size N m 10 4 . 
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All situations that lead to change are shown below: 


where a single-line arrow represents the state of a 
node belonging to the q- lobbies: arrows in the nu¬ 
merator represent states of the nodes belonging to 
the g-lobby chosen on first level and arrows in the 
denominator those on second level; voter’s state is 
marked by a double-line arrow. 

We should commend here, that the (/-voter model 
with the GLOBAL&AND rule and q = (/2 on the duplex 
clique is equivalent to the q -voter on the monoplex 
clique with q = q 1 = 2(/2 (i.e. a monoplex clique 
with a (/-lobby size twice as large as in the duplex 
case). This is visible in Fig. 0] where we compare 
the Monte Carlo simulations obtained for the duplex 
with those obtained for the monoplex topology. 

(ii) GLOBAL&OR - global independence and OR rule (see 
an example in Fig. [5]). 

Here, similarly as in the GLOBAL&AND rule, the status 
of independence is the same for both levels. With 
probability p the voter is independent and changes 
its state to the opposite one with probability 1/2. 
With probability 1 —p the voter behaves like a con¬ 
formist and its state is dependent on both q- lobbies. 
In contrast to AND dynamics, now the voter changes 
its state to the opposite one even when only one 
q -lobby is self-consistent and the second is not. In 
the situation when two (/-lobbies are homogeneous 
but not in agreement, i.e. one q -lobby supports the 
voter and second suggests to change its state, the 
voter becomes confused and stays in its old state: 
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(iii) LOCAL&AND - local independence and AND rule (see 
an example in Fig. [6]). 

In this case the independence is related to the 
level,i.e. we run dynamics separately on each level. 
It means that a voter is independent on the first 
level with probability p and with probability 1 — p 
behaves as a conformist — it is under the influence 
of the (/-lobby on this level. The same situation 
is on the second level, where regardless of the first 
level we choose if the voter behaves independently 
or conform the q -lobby on the second level. Finally 
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FIG. 5. The GLOBAL&OR rule: a voter is independent regard¬ 
less of the level with probability p and is subjected to the 
peer pressure with probability 1 — p if at least on one level 
the g-panel is self-consistent. In this example on level 1 the 
(/-lobby (agents in circles) is homogeneous and has the state 
opposite to the state of the voter (an agent in the square). 
Simultaneously, on level 2 the q -lobby is not self-consistent. 
Therefore the voter is not confused by two opposite (/-lobbies 
and is influenced by the first (/-lobby. 


we change the state of the voter only when both 
separated dynamics give in result the same state. 

We do not consider the LOCAL&OR rule because of the 
difficulty of such a concept related to social unreality and 
algorithmic ambiguity. 


IV. THE TIME EVOLUTION 


The aim of this section is to derive equations that de¬ 
scribe the time evolution of the system for each of three 
considered rules. Let us denote by N^(t) the number of 
voters in + 1-state (up-spins) at time t and by N±(t) the 
number of voters in —1-state (down-spins). The total 
number N of spins in a system does not change and we 
can define the concentration of up-spins at time t as: 


c(t) 


Ml 

N 


(7) 


Since all individuals keep the same state on both levels 
and we consider a duplex clique, we can simplify our 
analysis by considering concentration c(t) only on one 
level. However, we need to stress that the changes of the 
state of the node occur under the influence of both levels. 
In a single time step A t three scenarios are possible — 
the number of up-spins N^(t) will either: increase by 
1, decrease by 1 or remain constant. Simultaneously the 
concentration c(t) increases or decreases by or remains 
constant: 


7 +(c) = Pr {c(t + A t ) = c(t) + -1}, 


(8) 


FIG. 6. The LOCAL&AND rule: a voter is independent sep¬ 
arately on each level with probability p and subjected to 
peer pressure separately on each level with probability 1 — p. 
It means that on each level independently the voter (in the 
square) can be in one of three possible states: +1 with prob¬ 
ability p/2, —1 with probability p/2 and in a state suggested 
by the (/-lobby (in circles) with probability 1 — p. In this 
example, the (/-lobby is homogeneous on the first level and 
therefore it influences the voter. On level 2 the (/-lobby is not 
self-consistent and therefore there is no peer pressure. Finally, 
there are nine possible pairs that represent states on the first 
and on the second level - see the numbers in the corners of 
squares in the bottom line. The voter changes its state only 
if states, obtained independently on each level, are the same 
- see all possible final states of the voter in the bottom line. 

7“(c) = Pr {c(t + A t ) = c(t) - F} ; 

7 °(c) = Pr{c(7 + A t ) = c{t)) = 1 - 7 + (c) - 7 “(c). 

The time evolution of the average concentration is 
given by the rate equation: 

< c(t + At) >=< c(t) > +F [ 7 +(c) - 7 “(c)] , (9) 

where the exact formulas for probabilities 7 + (c) and 
7 ~(c) depend on the applied rule. In the following part 
of the paper, we use the abbreviated notation replacing: 
7 + (c) by 7 + , 7 - (c) by 7 “ and c(t) by c. Explicit forms 
of probabilities 7 +, 7 “ are the following: 

(i) GLOBAL&AND 

7 + = (1 - p)( 1 - c)c 2q + p( 1 - c)/ 2 , 

7 _ = (1 — p)c(l — c) 2q + pc/2. (10) 

The first component describes conformity, where 
the change of the state is possible only when two 
lobbies of size q each (i.e. the total number of agents 
is equal to 2 q) possess the opposite state than the 
state of the voter. The second component is respon¬ 
sible for the change due to the independence. 
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(ii) GLOBAL&OR 
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Since in the OR case agreement just in one lobby 
is needed, the sum in the above equations reflects 
all possible states in which all agents in the g-lobby 
on the one level possess the state opposite to the 
state of the voter and simultaneously the g-lobby 
on the second level is not homogeneous. GLOBAL&OR 
rule for q = 2 indicates change of the voter’s state 
if three or four of the four agents (two from each 
level) posses the same state and therefore this it is 
equivalent to the majority rule (35l-[38j. For q > 2 
majority rule is not enough since changes are pos¬ 
sible only when at least one lobby is homogeneous, 
e.g: 
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This fact is direct reason why in the following ex¬ 
ample there is no change of state: 
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(iii) LOCAL&AND 


7 + = (1 - p) 2 { 1 - c)c 2q + p( 1 - p){ 1 - c)c q + - ^ , 

2 

7“ = (1 - p) 2 c{l - c) 2q + p( 1 - p)c( 1 - c) q + ^. (14) 



FIG. 7. The average magnetization < m > as a function of 
the stochastic noise p for the GLOBAL&AND rule on the duplex 
clique. Monte Carlo results (empty symbols) were obtained 
for the system of size N — 10 4 and averaged over 10 3 sam¬ 
ples. The numerical solutions of Eq. © are marked with full 
symbols. 


V. RESULTS 


Solving analytically rate equation Eq. ([9]) in general 
(i.e. for arbitrary q) is a difficult task. However, it is 
easy to obtain a numerical solution by iterating Eq. ©• 
In such a way we can obtain the time evolution of the 
average concentration < c(t) >, as well as the stationary 
value < c >. The magnetization m(t) defined by Eq. ([2]) 
is directly related to the concentration c(t): 


m(t ) 


Ml z Ml 

N 


= 2 c{t) - 1 


(15) 


and therefore using the rate equation m we can easily 
find also the average magnetization. Independently, we 
can obtain results conducting Monte Carlo simulations 
and calculating the ensemble average of the magnetiza¬ 
tion defined by Eq. (|3]). 

Relations between the average magnetization in the 
stationary state < m > and the stochastic noise p, ob¬ 
tained by two methods (numerical solutions of analytical 
formulas and Monte Carlo simulations), are presented in 
Figs. [7J [9] and [51 It is seen that the agreement between 
the Monte Carlo results obtained for system size N = 10 4 
and the numerical solution of Eq. for the infinite sys¬ 
tem size (N M oo) is very satisfactory. 


Here we have three components: the first describes 
the situation when the voter behaves like a con¬ 
formist on both levels, the last one corresponds to 
the case where on both levels the voter is indepen¬ 
dent. The second term in Ea. (fTT|) is a mixed one - 
the voter behaves as a conformist ((1 — p)c q ) on one 
level and simultaneously it is independent on the 
second level (p/2). We multiple this middle term 
by 2 since this situation can appear in two configu¬ 
ration: conformist on the first level and independent 
on the second and vice-versa. 


VI. PHASE TRANSITIONS 

Due to the Landau theory, to describe any kind of 
phase transition we can introduce the quantity that mea¬ 
sure the degree of order (order parameter) [39|, [40[ . Al¬ 
though originally Landau theory was created to describe 
continuous phase transitions [39|, it occurred that the 
theory can be used also in the case of discontinuous phase 
transitions [io|. An order parameter, introduced to dis¬ 
tinguish between two phases, is equal to 1 in the com¬ 
pletely ordered state, decreases as a function of the devi- 














7 




FIG. 8. The average magnetization < m > as a function of the 
stochastic noise p for the GLOBAL&OR rule on the duplex clique. 
Monte Carlo results (empty symbols) were obtained for the 
system of size N = 10 4 and averaged over 10 3 samples. The 
numerical solutions of Eq. © are marked with full symbols. 



FIG. 9. The average magnetization < m > as a function of the 
stochastic noise p for the LOCAL&AND rule on the duplex clique. 
Monte Carlo results (empty symbols) were obtained for the 
system of size N = 10 4 and averaged over 10 3 samples. The 
numerical solutions of Eq. © are marked with full symbols. 

at ion from the order, and becomes zero in the disordered 
phase. Therefore for our system the natural choice of the 
order parameter is an average magnetization < m >. 

It is seen in Figs. 0 [9] and [8] that for all three rules 
GLOBAL&AND, GLOBAL&OR and LOCAL&AND, the system un¬ 
dergoes the phase transition. Below the transition point 
p = p c , the average magnetization < m > (order param¬ 
eter) is non-equal zero and above the transition point 
< m >= 0. For GLOBAL&AND (see full symbols in Fig. 0) 
the transition changes its type from continuous to dis¬ 
continuous at q = 3, which corresponds to q = 6 for the 
mo nop lex clique and thus agrees with results obtained 
in [13[. Analogously, for GLOBAL&OR (see full symbols in 
Fig. [8]) the transition changes its type also at q = 6. 
However, for LOCAL&AND (see full symbols in Fig. [9]), the 
transition changes its type from continuous to discontin¬ 
uous already at q = 5. Moreover, the transition point 
is much higher than for remaining two rules. The aim 
of this section is to determine the relation between the 
threshold value p c and parameter q for all three rules 
and understand deeper the nature of the observed phase 


FIG. 10. The relation between the critical value of the 
stochastic noise p c and parameter q for different variants of 
the model. 

transitions. 

We can obtain the the transition point p = p c directly 
from the Landau definition of an order parameter < m > 
using numerical solutions of equation m or Monte Carlo 
simulations. Because, as seen in Figs. 0 [HI and [9j the 
agreement between Monte Carlo and numerical station¬ 
ary solutions of Eq. m is very good, we determine p c 
from Eq. m, which is not only much faster but also more 
accurate method. Figure [TOl shows the relation between 
the critical value of noise p c and the size of the g-lobby. 

For the monoplex network the critical value of the noise 
for q < 5 has been derived analytically jl3-]: 

Pc(q) = --—-—r, (16) 

q — 1 -\- 2 q ~ x v ’ 

which gives p c ( 2) = p c { 3) = 1/3 and p c (l) = 0, i.e. there 
is no phase transition for the linear voter model (q = 1). 
However, for the multiplex structure the phase transition 
is observed even for q = 1, since we choose one neighbor 
from each level. Moreover, for q = 1 the critical value 
of noise p is equal for GLOBAL&AND and GLOBAL&OR rules. 
This is clear, because in this case, for both rules, the 
change of the voter’s state can happen only if the neigh¬ 
bors chosen from the first and the second level possess the 
same state, opposite to the voter’s state. For GLOBAL&AND 
the relation p c = p c (q) is a decreasing function for the 
whole range of g, which is also understandable, because 
this case corresponds to the q- voter on monolex with with 
twice the size of the q-lobby. 

From the point of view of the network structure, cer¬ 
tainly the most interesting, among considered rules, is 
the LOCAL&AND rule. For both GLOBAL rules multiplex 
could be in fact replaced by the monoplex network. In 
the case of GLOBAL&AND, as already mentioned, we could 
simply consider the g-voter on monoplex with doubled 
size of the g-lobby. The GLOBAL&OR is less trivial but still 
could be probably reformulated in terms of the g-voter 
model with the threshold on monoplex [20|. The case of 
LOCAL independence is not only less trivial, but also the 
most interesting from the social point of view. It should 
be remembered that conformity (and simultaneously in¬ 
dependence) is relative, i.e. individuals always conform 
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in respect to the particular social group and there are 
many factors that influence the level of conformity [El. l4ll— 
l43l ]. It means that the same individual may conform to 
one group and behave independently in respect to an¬ 
other. For example it has been shown, on the basis 
of various social experiments, that the level of confor¬ 
mity is much higher in the face-to-face condition than in 
computer-mediated communication gi Hi. Hence the 
idea of local independence is highly justified in modeling 
social systems. 


Therefore, we will now concentrate on LOCAL&AND rule 
and discuss the phase transition more thoroughly in this 
case. First of all let us notice that 

F = 7 +_ 7” (17) 

can be treated, analogously as in 0 , as an effective 
force - 7 + drives the system to the state ’spins up’, while 
7 “ to ’spins down’. Therefore, inserting explicit forms 
of 7 + , 7 _ from Eq. m we calculate also an effective 
potential: 


V = - J Fdc 

= p{ 1 - p) {<? [(1 - c)c 9_1 + c( 1 - c) 9_1 ] - [c 9 - (1 - c) 9 ]} 

+ (1 - p ) 2 {2 q [(1 - c)^- 1 + c( 1 - c) 29 ” 1 ] - [c 29 - (1 - c) 29 ] } 

- p 2 /2. (18) 


To find the critical value of the stochastic noise p c and 
the threshold value g, above which the transition becomes 
discontinuous, we could now use the Landaus approach, 
analogously as in jj3]. To do this we first rewrite the 
potential (fl8l) in terms of magnetization m, using rela¬ 
tion (usd, and then expand it into power series around 
m = 0. Unfortunately, the form of the potential is much 
more complex in this case than for the q -voter model on 
monoplex |l3|. Therefore all formulas are much longer 
and difficult to analyze. To understand the nature of the 
phase transition, it is much easier and more illustrative 
to draw potential V as a function of c for different values 
of p and q (see Figs. [Til and [T2| ) . 

For q <q — 5 (see Fig fTT]) the potential, given by Eq. 
m, behaves typically for the continuous phase tran¬ 
sition. Below the transition point the potential has 2 
minima that correspond to ordered states (i.e. 1/2 

and simultaneously m/0). With increasing p, minima 
are getting shallower and approaching each other. Even¬ 
tually they form a single minimum, that corresponds to 
the new disordered phase. For q >q the potential indi¬ 
cates discontinuous phase transition. For small values of 
noise p potential has 2 minima that correspond to ordered 
states (i.e. c / 1/2 and simultaneously m ^ 0). For 
larger p the third minimum (that corresponds to the new 
disordered phase) appears. Initially (for p E (Pi,^)) the 
minima that correspond to an ordered phase are deeper 
than the middle one , i.e. disordered state is metastable. 
For p = p 2 all three minima are equally deep, i.e. or¬ 
dered and disordered states are equally probable, which 
corresponds to the discontinuous phase transition. For 
p E (P 25 P 3 ) the potential has still 3 minima, but now two 
ordered states are metastable. Finally, for p E (p^l) 
the potential has only 1 minimum that corresponds to 
the disordered phase. Between spinodal lines i.e. for 
p E (Pi,Ps) one can expect hysteresis and indeed it was 


p=0.4 p=0.44 





FIG. 11. An effective potential V, given by Eq. m , as 
a function of concentration c for the LOCAL&AND rule and 
q — 4. For small values of noise p potential has 2 minima that 
correspond to ordered states (i.e. c / 1/2 and simultaneously 
m ^ 0). With increasing p, minima are getting shallower 
and approaching each other. Eventually they form a single 
minimum, that corresponds to the new disordered phase. This 
is a typical behavior for a continuous phase transition. 

found in Monte Carlo simulations (see Fig. [T3| ) . 


VII. CONCLUSIONS 

We have generalized the g-voter model with indepen¬ 
dence (noise) p to a duplex clique, i.e. a network that 
consists of two distinct levels (layers), each of which is 
represented by a complete graph (i.e. a clique) of size N. 
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p=0.3 



p=0.388 



p=0.38 



p=0.45 



FIG. 12 . An effective potential V , given by Eq. m , as a 
function of c for the LOCAL&AND rule and q — 5. For small 
values of noise p potential has 2 minima that correspond to 
ordered states (i.e. c ^ 1/2 and simultaneously m / 0 ). 
For larger p the third minimum (that corresponds to the new 
disordered phase) appears. Initially (for p £ (pi,P 2 )) the 
minima, that correspond to an ordered phase, are deeper than 
the middle one , i.e. disordered state is metastable. For 
p = p 2 all three minima are equally deep, i.e. ordered and 
disordered states are equally probable, which corresponds to 
the discontinuous phase transition. For p £ (P 2 iVV) the 
potential has still 3 minima, but now two ordered states are 
metastable. Finally, for p £ (^ 3 ,1) the potential has only 1 
minimum that corresponds to the disordered phase. Between 
spinodal lines i.e. for p £ (^ 1 ,^ 3 ) one can expect hysteresis 
and indeed it was found in Monte Carlo simulations (see Fig. 

m- 



FIG. 13. The average magnetization < m > as a function 
of the stochastic noise p for the LOCAL&AND rule and q — 5 
obtained from Monte Carlo simulations. Two different initial 
states were considered -’polarized’, i.e. ordered state with 
m — 1 and ’random’, i.e. disordered state with m — 0. As 
expected for discontinuous phase transitions (see Fig. fT2l) . 
hysteresis is observed. 


Levels represent two different communities (e.g. Face- 
book and school class), but are composed of exactly the 
same people - each node possesses a counterpart node 
in the second level. Such an assumption reflects the fact 
that we consider fully overlapping levels, being an ideal¬ 
istic scenario. We also assume that each node possesses 
the same state on each level, which means that the so¬ 
ciety consists of non-hypocritical individuals only. We 
have considered two criteria of level dependence - one 
related to the status of independence (GLOBAL vs. LO¬ 
CAL) and one related to peer pressure (AND vs. OR). 
The GLOBAL rule means that an agent is independent on 
both levels, but the LOCAL rule admits a situation where a 
person is independent in one clique but not in the other. 
Furthermore, the AND dynamics is more restrictive and a 
node changes its state only if both levels suggest changes, 
in the OR variant one level is enough to change the state 
of an individual. 

For all three considered rules (GLOBAL&AND, GL0BAL&0R 
and LOCAL&AND), the system undergoes a continuous 
order-disorder phase transition at p = p c (q ) for q < q 
and a discontinuous for q > q, where p c and q are rule- 
dependent. The GLOBAL&AND rule leads to a trivial result, 
identical with the monoplex case for a doubled value of 
q. For the GL0BAL&0R dynamics, p c is larger than for 
the monoplex network. However, q is identical with the 
monoplex case, i.e. q = 6 . In contrast to the other two 
rules, we find a qualitative change for the LOCAL&AND 
rule, as the phase transition becomes discontinuous for 
q = 5. The case of LOCAL independence is not only 
less trivial, but also more interesting and better justi¬ 
fied from the social point of view. In particular, it has 
been shown that the level of conformity during face-to- 
face communication is significantly higher than during 
computer-mediated communication such as the Internet 
0 , 0 . 

This suggests that the LOCAL&AND rule is the most suit¬ 
able for real social systems. Certainly it could be fur¬ 
ther developed by introducing different values of noise 
on each level. The simplistic duplex clique topology, as 
introduced in this paper, can be also modified to ob¬ 
tain a more general network. For instance, one could 
consider partially overlapping cliques, where some nodes 
possess no counter-node on the second level. Unfortu¬ 
nately, these modifications significantly complicate the 
model by introducing additional parameters and there¬ 
fore are beyond the scope of this paper. However, even 
considering such a simple model as here, we can observe 
that a multiplex network can introduce significant differ¬ 
ences in opinion dynamics. 
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